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GENERALIZED DYNAMIC ENGINE SIMULATION TECHNIQUES FOR THE DIGITAL COMPUTER 

James Sellers and Fred Teren 
NASA Lewie Research Center 
Cleveland, Ohio, USA 44135 


SUMMARY 

Recently advanced simulation techniques have been developed for the digital computer and used as the 
basis for development of a generalized dynamic engine simulation computer program, called DYNGEN, This 
computer program can analyze the steady state and dynamic performance of many kinds of aircraft gas turbine 
engines. Without changes to the basic program, DYNGEN can analyze one- or two-spool turbofan engines. 

The user must supply appropriate component performance maps and design-point information. 

Examples are presented to illustrate the capabilities of DYNGEN in the steady state and dynamic modes 
of operation. The analytical techniques used in DYNGEN are briefly discussed, and its accuracy is compared 
with a comparable simulation using the hybrid computer. The impact of DYNGEN and similar all-digital pro- 
grams on future engine simulation philosophy is also discussed. 

LIST OF SYMBOLS 
Symbols 

Ag exhaust nozzle area, m2 

A state matrix 

a coefficient 

Ei error variable 

f( ) function 
h enthalpy, J/kg 

J polar moment of inertia, kg-m2 

Ki control gain 

M matrix of aE^/aVj 

N rotor speed, rpm 

Ni high-pressure rotor speed, rpm 

N 2 low-pressure rotor speed, rpm 

P pressure, U/m2 

R gas constant, J/kg/K 

S Laplace transform variable 

T temperature, K 

C time,' sec 

u specific Internal energy, J/kg 

V component volume, m^ 

Vj independent variable 

<i mass flow rate, kg/sec 

y dependent variable 

Y specific heat ratio 

A incremental change 

c parameter in difference equation 

A eigenvalue of differential equation 

U eigenvalue of difference equation 

t time constant, sec 

-/- 
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«• 

state matrix 

Subscripts: 

accel 

acceleration schedule 

c 

compressor 

decel 

deceleration schedule 

dem 

demand 

F 

fuel flow 

i 

integer 

In 

into control volume 

j 

integer 

max 

maximum 

min 

minimum 

n 

integer 

out 

out of control volume 

r 

reference 

T 

turbine 

0 

base value 


INTRODUCTION 

The Intent of this paper Is to discuss new techniques in use at NASA Lewis Research Center for the 
simulation of turbojet and turbofan engine dynamics. An introductory discussion will be given on the rela- 
tive merits of analog, hybrid, and digital computers for use in dynamic engine simulation, but the body of 
the paper discusses a new digital computer program, called DYNGEN, which possesses significant advantages 
over traditional digital simulation methods . ! 

Computer programs which predict the performance of real and proposed aircraft engines have long been 
recognized aa indispensable tools for preliminary and detail design work. As engines and aircraft have 
grown more complex, analytical techniques have grown along with them to assist in optimizing engine con- 
figurations, comparing predicted performance with applicable criteria, exploring off-design performance, 
developing control modes and schedules, and providing timely data inputs to airframe designers. Recently, 
the prediction of engine dynamics has begun to play a significant role even in the preliminary design 
phases. For V/STOL aircraft, in which engine thrust provides lift and attitude control, engine response 
considerations- may have a decisive effect very early in design. Therefore, dynamic engine simulations may 
provide data for the most fundamental design decisions as well as fulfilling their traditional role of sup- 
porting control modes studies and other kinds of development experiments. 

The analog computer is one of the traditional tools of the dynamics and control specialist. Its chief 
advantage lies in the use of amplifiers which directly integrate the differential equations used to model 
the dynamic system. The analog computer is also a parallel-processing device, which means its components 
are all generating solutions simultaneously instead of one step at a time as in a digital computer. There- 
fore, the size of an analog problem has no effect on its execution time, and real-time solutions can some- 
times be obtained. However, experience has shown that large analog simulations can be unwieldy to operate 
because of hardware limitations and lengthy setup procedures. Some other problem areas in analog simula- 
tion are the difficulty of generating bivariant functions such as compressor and turbine maps, poor day- 
to-day repeatability of solutions, and the difficulty of transferring simulations to other users. Since 
analog computers cannot easily solve Implicit algebraic equations, simplified models often must be used to 
obtain explicit solutions for all variables (ref. 1). As a result, analog solutions will not agree per- 
fectly with the solutions generated by highly-detailed steady-state digital simulations. 

Modern hybrid computers have alleviated some of the problems noted above for analog computers. Specif- 
ically, the problem of bivariant function generation can be handled on the digital part of the hybrid com- 
puter, Also, the digital computer can be used to automate many of the setup and checkout procedures which 
confront the user of all-analog equipment. Although the digital part of a hybrid computer is slow compared 
with the analog part, real time solutions are still an attractive possibility for hybrid simulations. 

As digital computers have grown larger and faster their attractiveness for dynamic engine simulation 
has improved. Their chief advantages lie in their ability to solve large numbers of complex algebraic 
equations and the ease with which logic can be implemented. Bivariant functions can be simulated with data 
tables, and systems of implicit algebraic equations can be solved by iterative methods. Digital programming 
Languages have become standardized sufficiently to allow transfer of "engine decks" among users. Also, 
digital computers produce highly repeatable results. The main disadvantage of the digital computer, as 
far as dynamic simulations are concerned, is the need to choose an approximate method for solving the dif- 
ferential equations which model the system. Traditionally, this has led to numerical stability problems 
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and very long execution times for engine simulations. Furthermore, since most methodB for solving differ- 
ential equations require an explicit solution for all the derivatives (just as on the analog computer) , the 
analyst frequently must use equations which differ from those used in a purely steady-state program. As a 
consequence the digital dynamic engine simulation will often share a problem with Its analog counterpart: 
it won't agree with the steady-state deck. 

Clearly, traditional methods of dynamic engine simulation involve problems which are worth eliminating. 
The purpose of this paper is to describe an all-digital computer program called DYNGEN which successfully 
solves some of the problems noted above. Historically, DYNGEN is a derivative of the GENENG program 
(ref 8. 2 and 3) developed at NASA Lewie Research Center. GENENG, in tutu, is a derivative of the SMOTE 
program (refs. 4 and 5) developed by the U.S. Air Force Aero Propulsion Laboratory. SMOTE and GENENG are 
purely steady-state programs, but DYNGEN uses a recently developed method of solving differential equations 
which extends the capability of GENENG to include engine dynamics. As a result, DYNGEN Includes in a 
single program all the eteady-state capabilities of GENENG plus the added capability for dynamic calcula- 
tions. This eliminates the traditional discrepancy between the results of dynamic and steady-state pro- 
grams. DYNGEN is a generalized program: it enables the user to analyze one-, two-, or three-spool engines 

without reprogranming. Only component performance maps and certain design point information need be pro- 
vided. DYNGEN also allows a variable time step in computing the time -dependent solution of the differential 
equations used in the engine model. Thia feature shortens execution times for long transients where the 
user is concerned only with low-frequency dynamics. DYNGEN is written in Che FORTRAN IV language for the 
IBM 7094 computer. The basic version of DYNGEN requires about 32,000 words of storage. 

The description of DYNGEN will begin with an overview of the amount of detail included In the thermo- 
dynamic and component calculations. Next, a description of the procedure used to obtain steady-state op- 
erating points will be given. The discussion will then proceed to explain how the solution of differential 
equations can he made a natural extension of the steady-state techniques by using a modified Euler method 
for solving differential equations. Appendix A is included to give mathematical details of the Nevton- 
Raphson iterative technique, which DYNGEN uses to obtain bnth steady-state and dynamic operating points, 
and the modified Euler method for solving differential equations. 

The latter portion of the paper is devoted to user-oriented subjects. A few examples are given to 
show the variety of engines that can be simulated without reprogramming, and the possible options for spe- 
cifying off-design points are described. Finally, some examples of transient operation are presented to 
show how DYNGEN operates in connection with user-supplied control system subroutines. 


ANALYTICAL TECHNIQUES 

Thermodynamic and Component Calculations 

Since DYNGEN is a modified version of a steady-state program, it contains details which are not always 
found In purely dynamic simulations . A discussion of these details will be given to assist the reader In 
visualizing the level of sophistication In the analytical model without actually presenting specific equa- 
tions . 

Engine components are represented in the usual way by performance maps. A typical, fan or compressor 
map is shown in figure 1: pressure ratio and efficiency are plotted as functions of corrected flow and 
corrected speed. A typical combustor map is shown in figure 2: efficiency is plotted as a function of 
temperature rise across the combustor and entry pressure. The form of the turbine maps is one which has 
become Common in steady-state engine simulations. A work parameter, Ah/T, and efficiency are plotted as 
functions of corrected speed and flow parameter (iV?/F, aa shown in figure 3. Afterburner efficiency is 
calculated using three separate functions, as shown in figure 4. The user specifies a design-point after- 
burner efficiency, which Is then adjusted to account for changes in fuel-air ratio, afterburner Inlet Mach 
number, and afterburner inlet total pressure. 

DYNGEN has been programmed to provide automatic map scaling. This feature 1 b UBeful for preliminary 
design work since it means the user need only supply maps which he thinks would resemble the true maps of 
the engine he Is simulating. Xn addition to supplying maps, the user must also specify the design operat- 
ing point for each of the maps. The program will then linearly scale the map inputs and outputs to make 
their values compatible with the design pressure ratios, flows, etc., which have also been specified at the 
design point. The scale factors computed at the design point are saved and applied to the map inputs and 
outputs for all off-design cases. 

In all thermodynamic calculations, gas properties are calculated as a function of temperature and fuel- 
air ratio. The temperature rise across a compressor is calculated without resorting to assumptions about 
the average gas properties across the component. Instead, a three-step procedure is used: firBt, the 

isentropic enthalpy rise is computed from knowledge of the pressure rise; second, the isentropic enthalpy 
rise is corrected by the efficiency; and third, the actual temperature Is calculated from the actual en- 
thalpy and pressure. A similar procedure is used to compute the temperature drop across a turbine. 

Turbine cooling bleed is accounted for by mixing the bleed air downstream of the turbine; the bleed 
itself is assumed to do no work. In its basic form, DYNGEN contains no provision for compressor variable 
geometry or Interstage bleeds. This problem was left to the user to account for by making appropriate ad- 
justments to the component maps. Where the core and duct streams come together in a mixed-flow turbofan, 
a static-pressure balance calculation is made. The user also has the option of specifying separate noz- 
zles for the core and duct if he wishes to simulate an engine with no flow mixing. Finally, in line with 
current NASA policy, DYNGEN was written to accept either English or SI units. This is accomplished by 
changing physical constants within the program rather than converting only the input and output. 
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Steady-State Balancing Technique 

An example case will now be presented to assist the reader in understanding how DYNGEN calculates en- 
gine steady-state operating points. For simplicity a turbojet engine will be used in Che example, but 
similar methods are used for more complicated configurations. Figure 5 shows a turbojet engine with its 
major components labelled. Pressures (F) , temperatures (T) , and flows (w) are also labelled with appropri- 
ate station numbers. The example will illustrate how the calculation of Variables proceeds through the 
engine. DYNGEN is written so that the user can select off-design points by specifying speed (N) , turbine . 
inlet temperature (T4), or fuel flow (up). In this example, fuel flow is assumed to be the specified var- 
iable. First, an inlet calculation 13 made to determine F2 and T2 from the free-stream values of pres- 
sure, temperature, and Mach number. To calculate £ c , P3, and 1 3 from the compressor map (fig. 1 ) and 
thermodynamic relations, guesses must be made for the values of speed (N) and pressure ratio (P3/P2). 

Once (j c , P3, and T3 are obtained, the combustor calculations for 004, P4 , and T4 can be made using 
thermodynamic relations, the combustor map (fig. 2), and specified values for fuel flow (ciip) and compressor 
bleed flow. To calculate turbine variables another guess is made, this time for the value of turbine flow 
function (i>4 . Then, from the known value of (N/ yTT) , the turbine map (fig. 3 ) is used to calculate 
turbine work (Ah) and efficiency. P7 and Ty are then calculated using the thermodynamic relations. 
Finally, the compressible-flow relations are used to calculate nozzle pressure (P7) from log, T7, and spe- 
cified values for P 0 and nozzle area. 

The reader may have noticed that the above calculation procedure is redundant; that is, certain vari- 
ables can be calculated In more than one way. This fact Is used to generate error variables which must 
equal zero to yield a consistent solution of the simulation equations. In writing a program such as 
DYNGEN, the analyst has great freedom in choosing what error variables to use. This discussion simply 
points out the choices which were made by the authors of SMOTE; experience has shown that these were good 
choices for most engine configurations, and the same error variables were retained in SMOTE' s descendants, 
GENENG and DYNGEN. 

In the previous discussion it was stated that guesses were made for rotor speed (N) , compressor pres- 
sure ratio (P3/P2) , and turbine flow function (i&4 VT4/P4) . From the first two guesses (and other variables) 
one may calculate the power absorbed by the compressor (tO(.Ah c ) . From che turbine flow function (and other 
variables) one may calculate the power supplied by the turbine (wy-Ahy) . For steady-state Operation the 
power supplied must equal the power absorbed. Therefore, the difference (ii c Ah c - di^Ahi) may be used for an 
error variable. Similarly, one can calculate a value for turbine flow function ((114 ' based only on 
the first two guesses, but for a consistent solution the calculated value must equal the guessed value. 
Hence, the difference (014 - (ill 4 ' can be used as a second error variable. Finally, from the 
ccmpressible-f low equations we know that the variable Py is specified by the variables tig, Ty, P Q , and 
nozzle area. This value for Py must equal the value F j which is calculated from thermodynamic relations 
at the turbine exit and from pressure losses In the duct between turbine and nozzle. Therefore, the third 
error variable is (F7 - p )) . 

Once three variables hove been guessed and three errors have beer) specified, the analyst can use an 
iterative method to obtain a consistent solution to the simulation equations. SMOTE, GENENG, and DYNGEN all 
use the Newton-Raphson technique of iteration. The details of this method are given In Appendix A. In 
figure 6, a simplified flow chart shows how the Newton-Raphson method is used in connection with the en- 
gine calculations discussed in the preceding example. Although more complicated engines will require more 
guesses and more error variables in the iterative procedure, the analysis will be quite similar to the one 
described in the example. 


Simulation Differential Equations 

So far the discussion has been devoted to the methods which DYNGEN uses to obtain steady-state operat- 
ing points. Now the method of implementing and solving time-dependent differential equations will be dis- 
cussed. DYNGEN uses a modified Euler method of solving differential equations. This method is derived from 
a numerical analysis viewpoint in Appendix A. Appendix A also discusses the numerical stability of the modi- 
fied Euler method and shows that it does not require extremely small time steps to obtain a stable solution. 
This advantage is important in engine simulations because in the past it has often been necessary to select 
integration time steps small enough to guarantee stability for high-frequency dynamics typical of mass and 
energy storage in unsteady flow. This can result in very long execution times even though the simulation 
user may only be interested In low-frequency dynamics. With the modified Euler method, the user can select 
longer time steps without worrying about numerical stability. The main disadvantage of the modified Euler 
method is that an iterative solution is required for the difference equations which approximate the solution 
to the differential equations. However, this fact turns out to be useful in DYNGEN since it means that the 
analyst no longer has to solve explicitly for derivatives. They may be embedded anywhere in an overall set 
of simultaneous algebraic equations which are. to be solved by an iterative method such as Newton-Raphson. 

The following discussion shows how this advantage was employed in modifying a steady-state simulation, GEN- 
ENG, to form a dynamic simulation, DYNGEN. Figure 7 ehows the three kinds of equations which were modified 
to include dynamic terms: the power balance, continuity, and energy equations. The steady-state power bal- 

ance equation simply Implies that the power output of a turbine must equal the power absorbed by a fan or 
compressor. By adding a rotor acceleration term, che equation can be used to model engine dynamics: any ex- 

cess power provided by the turbine will go into rotor acceleration. If the time derivative is arbitrarily set 
equal to zero, the dynamic equation becomes the steady-state equation. Similar considerations also hold for 
the continuity equation, DYNGEN treats unsteady flow dynamics in a way which has become traditional for en- 
gine simulation: a control volume is associated with each component, and pressure, temperature, and density 

are assumed constant throughout the control volume. At steady state the flow into the volume must equal 
the flow out, but for unsteady flow mass can be stored in the volume at a rate proportional to the time 
derivative of pressure, dP/dt. If dP/dt is zero, the continuity equation reverts to its steady-state 
form. The control volume approximation is also used for the energy equation. At steady-state the rate of 
energy into the volume must equal the rate Out, but, in unsteady flow, energy storage is accounted for by 
two terms: one reflecting the rate of change of specific internal energy, du/dt, and another reflecting 
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energy storage caused by mass storage. 


DYNGEN was formed from GENENG by mo 
ever they occurred in GENENG. In GENENG 
ror variable: 


difylng the power balance, continuity, and energy equations where- 
, the steady-state power balance equation was used to form an er- 


^1 - — oyfQfi'p 

In DYNGEN, the same error is formed with the dynamic term added: 


E 1 = «c Ah c ~ «T flh T + JN ~ 


“ J '” ISy ™« *ss.cl»ted »£«, each 

if a, w a. fi» ;L:s; “ss ^T 1 " 1 by the p “ —pi- 

" th * c °* pr ““ r th “ the n ™ 

.i . . V3 dP3 

c c YRT3 dt 


h 3 


“c h 3 - (u c 


w 0 )u-j 


P3V3 du3 

RT3 dt 


The derivatives are calculated by the simplest possible approximation: 


i2. . yi - yi-i 

dt At 


approximation ^ “ad equate provided thl^th^timfs tep^A r for the P^vions time step. This 
time constant the user wants to observe Fer ule’ if ’ 'Ll ^-teath the magnitude of the smallest 

1.0-second time constant, he should use*a At no greater tLTIo’SSnJ? 4ynamlCS W±th 3 

ation^c^^^lS^ 6 ™^: Sh / 11 n0t reqUlre aUy Chan8e t0 the basic ite - 

capability was simply extended to include dynamics. * lblltty ° r e eneral ity of the program was lost; its 


DYNGEN CAPABILITIES 
Engine Configurations 

ihe r2L^;^T I £LS B fl ss B 22'e!: l\SK2:d a V e ir qoe F T d in DYNGEN > or hou “'«*»• 

a few examples will be given to show the varietv of Wh ? C Ca ” be done with ltl F±rs t. 

These examples are meant only to suggest the rakge of cot-ions th ^ be ana ^ 2ed without reprogramming. 
“°" 8 ’ , c W nt« lf.t „f £' 'too «n|loe =o«fl|o„- 

Chree-epoo2 J 6 three-streJ C ’eogioe! li rhe e orUlin nn" i ^ , the u 48lC ver8 * 0[l of DYNGEN: a 

to study blown-flap or elector-wing orooulslon <n V «t f ™?, 3 ® uiatin 8 this kind of engine was the need 
the "wing duct," and its aL doL not mix wJth III f 0L aircraft. Hence, the third duct is called 

be eliminated to simulate a three-spool two-stream^nod nl ^ UCt ® traani ®- If deslred » the wing duct can 
another step down in complexity: a two-SDool hn-> tf>dff SU &S *55 R ° ll8 ~ Ro Y cs ^-211. Figure 9 shows 

by separate component maps. All of the turbofan ■ an en S ne - The fan and booster may be represented 
rather than separate core ducf noSle^s sho^^n h « by D¥NGEN C “ 1186 a option 

also^avanabie options. Pinally, fi g nre 10 shows the 


Steady-State Capabilities 

predecessor, GENENG. . For^Lr^gL^configuratiorthlrthru^r^i 1 ^ 163 U ' llCh DYWGEN inherited from its 
must be run. At the design point the user must specify certain cycle market 468180 P ° lrlt C3Se 
perature, component efficiencies, Mach numbers at various " cycLe parameters such as turbine inlet tem- 
must be provided, and logical controls must be set to establish how etc - Con P°n a nt maps 

core and duct streams mix, whether it is a turboiet or turhnfan h , any s P ooFs c he engine has, whether its 
ponent maps will be automatically scaled to conform with the^M." 4 S ° ° n ‘ ^ tnenCloned e «lier, the corn- 

design point. After the design point case is run DYNGfn u/ 1 param f te * s specified by the user at the 
burner fuel flow, thrust area! a l vario^ a^tion^ fZTrl ^ ^ »**. — after- 

First, there are four basiroperating^odes for siecif5ing f off-desi Stead r State P ° intS (refS ‘ 3 and 4) • 
perature, constant main fuel flow, constant fan j V 8 ff de8±gn polnts- constant turbine Inlet tem- 
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ating parameters such as altitude, Mach number, inlet recovery, bleed or power extraction, and core, duct, 
or wing nozzle areas. If afterburning Is to be investigated, the user may specify either afterburner fuel 
fiow or temperature.. Nozzle area will be recalculated to maintain the same core engine operating point 
specified for some previous, nonafterburning case. These options also hold for duct burning. If the user 
wishes to investigate the effect of a converging-diverging nozzle on thrust performance, he may hold exit 
area constant or have it automatically recalculated for fully expanded flow. These options make DYNGEN a 
very useful tool for investigating steady-state off-design performance. 


Dynamic Calculations 

Thc f?v f r Exam P le % wil! now be given to show some of DYNCEN's capabilities for simulating engine dynamics. 

at exampie will show an open-loop engine response to a step in fuel flow, and the remaining two will 
demonstrate how DYNGEN can be used In connection with user-supplied engine control subroutines. The input 

cent th“t nt % f ° r r 8lne dynam1cs are identical t0 th « requirements for steady-state operation ex- 
cept that a few additional constants most be supplied at the design point: rotor moments of inertia, rotor 

design speeds in revolutions per minute, and component volumes. A design point case must be run, just as 
in steady-state operation, and if the user wants to specify some off-design point as the initial condition 
for a transient, he may do so using the steady-state options discussed previously. Finally, the user must 

8 timE ]f £ep £or “ansient solution, and he must set a program index indicating that a transient 
run. e transient disturbance input is supplied by a user-written subroutine. Any of the off- 
design input parameters which are available for steady-state operation, such as fuel flow, compressor bleed, 
inlet recovery, etc., may also be used as transient inputs. Furthermore, the user may generate/transient 
inputs appropriate to the control system he Is simulating, for example, a change in power lever angle. 

^ fl ^ St ex f npl f sbowa the response of a three-spooL, three-stream engine (like the one shown in fig. 8) 
to an open loop step in fuel flow. Figure 11 shows time histories of fan, middle spool, and core speeds. 

Also shown is the response of turbine inlet temperature. All variables are presented as percentages of 
their design values. Apart from showing DYNGEN 1 s capability to simulate a three-spool engine, figure 11 alsn 

eauaHoL atE5 Tb the ° f ^ tlTOe 8tepa in ChE solution of "he ^muSi™ ° 

equations. The results are shown for two time steps: 0.01 second and 0.10 second. Close examination 

shows some small differences between the two solutions, but they are substantially identical. There is a 
STin’ h ° We f r : ln “"P uter exec ution time to run the three-second transient shown in figure 11. 

Sme wls lS*T time step, execution time was 1.4 minutes; using the 0.01-second time step, execution 

rx' 3 mlnutes - ™ ls exam P le demonstrates one of the main advantages of a modified Euler solution 
method, the user may select the time step to show the frequency range of interest. If low- frequency ef- : 
8UCh 88 r0t ° r f dynamic*, are the subject of interest, a time step of 0.10 second may be adequate. If 
^ r ? h Hn f n ^ Cy , e ! feC£S ' SUCh 35 ten! P era ture and pressure dynamics, are to be observed, then a smaller time 
teresS 1 needed ‘ In any Case ’ e *ecution times can be held Co a minimum compatible with the user's in- 


ne ^ examp£e ehows f large throttle transient for a two-spool turbofan similar to the one shown in 
. Th± a gdne was simulated along with the speed control system shown In figure 12. The primary 
: t f r;S system is demand speed N 2idem which is set by the pilot's throttle lever. The only 

control L ironnrMo t ^ ! ° W ’ T ° h 5°*® C ° the C00 > bU8t0r - During small throttle transients the 

control is proport ional-plus-integral on speed error, but for large transients the control is closed-loop 
on the acceleration fuel flow schedule. Acceleration fuel flow is computed from compressor speed Ni com- 
pressor exit pressure, P 3 , and compressor inlet temperature T 2 ,i. This moderately complex control system 
was simulated using subroutines compatible with DYNCEN's modified Euler solution method! A throttle step 
from 50 percent thrust to 100 percent thrust was applied to the simulation, and the results are shoL 11 
figure 13. Time histories of thrust and turbine inlet temperature are shown with the variables expressed 

thos/f!^ Th°L thelr dSSiS T V f 1Ue ‘ ma f±gUre alS ° ^sents a comparison of DYNGEN ' s results wUh 
3 hybrld c °“P ute r simulation of the same engine. In figure 13, the continuous lines are the 
w computer solution and the discrete points are DYNCEN’s solution. The hybrid computer model is quite 

two%imJlltion 9 6 diff Ut T u £ ! differencea i" the simulation equations, the steady-scate results of the 
two simulations differ by about 3 percent. The differences in the dynamic solutions are of the same order 
The comparison shown in figure 13, though not perfect, tends to confirm the validity of DYNGEN's method of 
solving the differential equations used in modeling the engine and control system. Even though a fai!lv 

- —■ 

the °L D ™^' 9 flexlb / lity lnvolves a Single-spool, afterburning turbojet similar to 

the one shown in figure 10. This type of engine requires exhaust nozzle and main fuel control subsystems 
as shown in figure 14. The main fuel control is a simple proportional control on speed error with accel- 
eration and deceleration fuel flow limiting. The main input Is demand speed, H ’ which is set bv the 

L e ' ^ * cc f leratlon schedule is the usual (m F /P 3 ) accel as a function of N and T, and 
the deceleration schedule is obtained simply by taking one-third of the acceleration schedule. The'nozzle 
control is used only in the afterburning mode of operation. Its purpose is to null out any change in com- 
™rT re r I 2 Jl 2) ’ WhiCh mi8ht ° CCUr When 8°ing from non-afterburning to afterburning ope^- 

sir! , raSo S e ^o!! COn,P 7 ptoportional - plus ~ in ^g™l control of nozzle area, A 8 , in response L Jres- 

This control system was simulated in connection with a turbojet engine, and a throttle slam from idle 
to full afterburning was applied. The results are shown in figure 15. Time histories ^ rotational speed 
° fuel flow, afterburner fuel flow, nozzle area, and thrust are shown. All variables are presented^ ’ 
percentages of their design value. To simulate a throttle slam, afterburner fuel flow was ramped f!om zero 
to its -Ximum value in two seconds, beginning as soon as rotor speed reached 100 percent. This example 
shows that DYNGEN can be used euccessfully to simulate the dynamics of an afterburning engine Furthermore 
it demonstrates that DYNGEN is not limited to small-perturbation problems. It has all the capabilities of ’ 

P ZtTTT TS !° C si T lat£n S gr ° SS but it Is faster than most traditional simulations tL 

five second transient shown in this example required about two minutes of computer execution time. 
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CONCLUSIONS 

A generalized digital computer program for simulating the steady-state and dynamic performance of 
turbojet and turbofan engines has been described and discussed. This computer program, called DYNGEN, 
posa eases significant advantages over traditional methods of digital engine simulation. Specifically, it 
eliminates the need to operate two separate computer programs to obtain steady-state and dynamic results. 

It uses a modified Euler method for solving differential equations which enables the user to specify a 
solution time step compatible with the frequency range of interest. This saves computer execution time 
when long transients are to be run. Finally, DYNGEN can simulate a wide variety of engine types without 
reprogramming. This saves money and man-hours when new engines are to be simulated. Wien real-time engine 
simulations are required the analyst must turn to analog or hybrid methods to achieve his goals. However, 
owing to the new digital simulation techniques used in DYNGEN, all-digital engine simulations are now 
capable of accomplishing nearly all the user’s tasks with convenience and flexibility. 
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APPENDIX A 


ANALYTICAL TECHNIQUES 
Steady-State Balancing Technique 

The following discussion explains the Iterative method which DYNGEN and its predecessor GENENG use to 
calculate steady-state operating points. As noted earlier, the calculation of a steady-state operating 
point requires solution of a system of nonlinear equations, corresponding to various engine matching con- 
straints such as rotational speeds, air flows, compressor and turbine work functions and nozzle flow func- 
tions. In order to satisfy these constraints, there are available an equal number of engine parameters 
which may be varied, such as compressor and turbine pressure ratios and flow functions. The specific num- 
ber of engine parameters (independent variables) to be varied and engine matching constraints (dependent 
variables) to be satisfied depends on the type of engine configuration being studied, and varies from three 
for a single-apool turbojet engine to nine for a three-spool engine. The computer program (DYNGEN) 
searches for the values of the engine parameters which result in the engine matching constraints being 
sat isfied. 


If the independent variables are denoted by Vj and the dependent variables by E^, the matching equa- 
tions can be written as 


Ei(Vj) = 0 


i = 1,2,. ...n 
j = 1,2,. ...n 


Thia is a set of nonlinear equations, which must be satisfied for a steady-state solution. The procedure 
used to satisfy these equations is the multi -variable Newton-Raphson method (ref. 7). With this method, 
changes in E are assumed to be related to changes in V by first-order, finite-difference equations: 


AE - MAV 


where AV and AE are n-vectors denoting changes In V and E from some reference condition, and M is 
an nxn matrix of partial derivations of E with respect to' V: 


«ij 


3Ei 

avj 


The matrix M is obtained by calculating a reference case and n Independent perturbed cases, such that 
only the J-th variable Vj is perturbed from its reference value on the J-th case. Then for the J-th 
case, 

AEh 

* AVJ 1 - n 

Once the matrix M Is obtained, the reference case is improved by using 

V = V r - M~lE r 

If the system of equations were linear, this process would lead to convergence in one iteration. In prac- 
tice, nonlinearities in the system prevent immediate convergence. In this case, the new V and E are 
taken to be the reference values, and a new matrix is generated. If the system is not too nonlinear, and 
initial guess for V are reasonably accurate, convergence is achieved in several iterations. 


Dynamic Equations 

Once an initial steady-state solution has been obtained, a time-varying solution may be generated. 
This requires the solution of a set of differential equations which model the system. The specific equa- 
tions which are used to model the engine were discussed in the text. In this, section, the procedure used 
to solve the differential equations in DYNGEN will be discussed. 

Consider first the differential equation 


“ fCy,t) (1) 

In order to obtain a numerical solution using a digital computer, this differential equation must be re- 
placed by a difference equation, in such a way that the solution of the difference equation is In some 
sense close to that of the differential equation. There are many ways in which this can be done, as dis- 
cussed. for example, in reference 7. A common method is to use a difference equation of the form 

Yj+1 “ Yj + At[ef(yj,tj) + (1 - E )f(y ;)+1 ,tj +1 )] (2) 

where 

>’j “ y(t 0 + jAt) 
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0 < e < 


1 


The bracketed quantity In (2) represents a weighted average of the derivative f(y,t) over the integration 
interval [tj ,t j+]_] . For e = 1, equation (2) becomes 

yj+1 " >'j + Acf ^ y j' c j^ 

Equation (3) ia known as Euler's method, and allows explicit calculation of yj+i as a function of the pre- 
vious values yj and t j . On the other hand, for e f 1, equation (2) is theroodified Euler method, and 
in general cannot be solved explicitly for yj+i, because of the dependence of f -on >’j+l which appears 
on the right hand side of the equation. In this case, some form of iteration must be u9ed at each integra- 
tion step to solve for yj+i- 


From the standpoint of simplicity of the integration formula, use of (3) is clearly preferable to the 
use of (2). However, there are two other important considerations: accuracy and stability. As discussed 

in the literature (e.g., ref. 7), use of (2) can lead to greater integration accuracy. Even more important 
for the dynamic engine simulation problem is the stability consideration. 


To illustrate the stability consideration, consider the linear differential equation 


dt 


ay 


(4) 


For this equation, (2) becomes 


yj+1 “ yj + ait [eyj + (1 - e)yj+i] 
which can be solved for yj+i to give 

( 1 + aeAt \ 

y J +1 “ + aF.it - aAt) y j 

the general solution for yj can be written 

7 $ = tJy 0 

where 


(5) 


( 6 ) 


(7) 


1 + as At 
1 + aeAt - aAt 


( 8 ) 


the original differential equation (4) is stable for a < 0; the difference equation solution (7) is stable 
for jrj <1. From (8), the requirements for stability of (7) can be established in terms of the require- 
ments on integration step size. At. Solving (8) for At yields 


At 


1 - r 

a(er - r - c) 


(9) 


The upper and lower bounds for At 


are obtained by setting r " ±1 
2 .1 


At < 


a(l - 2c)* 


At - unconstrained, c < w 


in (9) . This results in 


(10a) 

<10b) 


In particular, for the Euler method (e - 1), the step size must be less than (-2/a) in order to avoid 
numerically-induced instability, while for e < 1/2, the numerical method leads to a stable solution for 
any value of integration step size. 

The above results are readily generalized to a system of linear differential equations. Consider the 
system of equations 


^ - Ay (11) 

where y Is an n-vector and A i8 the tixn system matrix* Use of the numerical algorithm in. (2) results 
in 


y-) + l ■= yj + AAC[cyj + (1 - Oyj +1 ] 

which has the general solution 


( 12 ) 



( 13 ) 


t = (I + AtAt - AAt)“l(I + AeAt) 



where 
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As shown in reference 8, (11) is stable If and only if the eigenvalues of A all have negative real parts 
the difference equation solution (13) is stable if and only if all the eigenvalues of <j> have magnitude 
less than unity. 

It will now be proved that if A is an eigenvalue of A, then 




1 + AeAt 
1 + AeAt - lit 


(14) 


is an eigenvalue of t. 


Proof 


Let A is an eigenvalue of A, Then 


If p is an eigenvalue of . t>, then 


(A - XI | a 0 


But 


|* - pi | » o 


But from (14) , 


| * - pi)' - f (I + AeAt - AAt)-l(I + At At) - pl| 

_ I (I + AsAt) - y(I + AeAt - AAt) | 

[i + AtAt - Mt| 

„ 1(1 - n) (1 + AtAt) + pAAtl 
| I + AtAt - AAt | 


1 - P 


-XAt 

1 + XeAt - XAt 


so that 


U _ t i = I -XAt (I + AtAt) + (1 + XeAOAtA 
(1 + AsAt — XAt) 1 1 + AeAt - AAt 


which completes the proof. 


At| A - Al| 

<1 + XeAt - XAt) 1 1 + AeAt - AAt) ” 0 


similarity of (14) and (8), together with the requirement that all eigenvalues p have magnitude 
less than unity, allow the conclusion, similar to (10), that * 


{ (1 - 2t) ’ 


(15a) 


At - unconstrained, e < — 


(15b) 


the'ste^size'ls'restricted'by ° f A ** Io P^tiCular, for the Euler method, 


w 

in order to avoid numerical Instability. 

r ! SUl n S "* Va < ld ^ f ? r * Unear s y stem * ™ such general proofs are available for 
ear y. terns. However, in an intuitive sense. It seems reasonable that equation (16) is applicable 
to nonlinear systems if the matrix A and eigenvalues X are interpreted as -average'' values over an in- 
tegration step, and the system of equations (11) is not too nonlinear. 

, . of . equa J io " ( 4 16) ’ Particularly for the dynamic engine simulation problem, is the fol- 

lowing. The dynamic engine simulation generally contains a mix of high and low frequencies. The hleh fre- 

n£ss C and 171 ^ ^ 1 “P ed " vo ^ UIne representation of component dynamics, which includes the storage of 

mass and energy. The low frequencies result, for example, from rotor dynamics, and the slow motion of ex- 
“""“If “* “«” 1 Wjic. Frequently , the al ™l.tlu» S.r 1. lutureatad £ S-L™.", 
^ as overall engine spool-up time, and is not concerned with high frequency effects. - Typical 
transients arc of five— to ten— seconds in duration r 

If the simulation uses Euler's method, the integration step size la restricted by the highest frequenc 

? VET’ CV ^ th ° Ugh t le USer iB n0t ^terested in high frequency information. In thi i TtlT 

size of 10 4 seconds, or smaller, ls f requefltly requlred . the other hand5 tf an iapUcit \Z d lkti P 
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Euler) technique is used (e < 1/2), the step size is not restricted. It can be chosen to suit the desired 
frequency content of the output, which typically allows a step size of 0.1 seconds or larger. 


ITERATIVE SOLUTION PROCEDURE 

A problem which exists with the use of implicit methods, as noted earlier, is that for nonlinear dif- 
ferential equations, some iterative scheme is required to solve for the values of y, + ^ at each integra- 
tion step. The differential equations corresponding to the dynamic model of the engine may be written as 

ft = f(y) d7) 

where y and f are vectors. The state vector y represents pressures, temperatures and rotor speeds. 
The dimension of y (and f) depends on the type of engine configuration being studied. Nine state vari- 
a lea are required for a single-spool turbojet engine, and a greater number for more complex engines. The 
number of state variables required for a dynamic solution always exceeds the number of steady-state itera- 
tion variables required for the Newton-Raphson iteration discussed earlier. 

The difference-equation representation used in DYNCEN utilizes e = 0, so that (17) becomes 

Tj+1 = Yj + £tf(yj +1 ) (18) 

The discussion of the sample configuration in the main body of the report showed how the dynamic equations 
are incorporated into the structure of the steady-state solution. The steady-state continuity, energy and 
power equations are modified to be dynamic equations. The resulting dynamic equations are then Included 
either as error equations, or used to calculate flows and enthalpies at various stations throughout the 
engine. & 
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Figure 1. - Fan or compressor map. 
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FUEL-AIR RATIO 


(a) GENERALIZED AFTERBURNER COMBUSTION 
EFFICIENCY AS FUNCTION OF FUEL-AIR RATIO. 
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Figure 4. - Afterburner maps. 
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Figure 5. - Steady -state engine calculations for a turbojet. 
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Figure 6. - Flow chart of balancing technique. 


POWER BALANCE: 

STEADY STATE 

*T Ah T ■ \ Ah c 
DYNAMIC 

dN 

w T Ah T “W c Ah c + JN-^ 

CONTINUITY: 

STEADY STATE 

“'OUT = *IN 
DYNAMIC 

V dP 

w OUT " W IN * yRT dT 

ENERGY: 

STEADY STATE 

W OUT h OUT B % h IN 
DYNAMIC 

A OUT h OUT = ^IN h IN " ( *|N - *0UT ,u 

PVdu 
'RT dt 

Figure 7. - Simulation dynamics. 
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Figure 10. - One- and two-spool turbojets. 
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Figure 11. - Three-spool turbofan; response to fuel flow step. 
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Figure 12. - Two-spool turbofan speed control. 
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NOZZLE CONTROL 

Figure 14. - Afterburning turbojet control system. 
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Figure 15. - Afterburning turbojet response to throttle 
slam. 


NASA-Lewis 


